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Abstract 

Recent experimental results on slow light heighten interest in nonlinear Maxwell 
theories. We obtain Galilei covariant equations for electromagnetism by allowing 
special nonlinearities in the constitutive equations only, keeping Maxwell's equations 
unchanged. Combining these with linear or nonlinear Schrodinger equations, e.g. 
as proposed by Doebner and Goldin, yields a Galilean quantum electrodynamics. 
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1 Introduction 

In this letter we propose nonlinear constitutive equations restricting the symmetry of 
Maxwell's equations to Galilean symmetry. We also show that these constitutive equa- 
tions arise as the formal nonrelativistic limit (taking the speed of light c — > oo) of the 
nonlinear relativistic theory: no additional hypotheses about relative field strengths are 
needed. Admitting such nonlinear constitutive equations leads to possibilities that are 
entirely new. Maxwell's equations in our approach stay unchanged; it is the choice of 
the constitutive equations alone that makes the difference between relativistic and non- 
relativistic theories. We stress that the resulting Galilean electrodynamics is essentially 
nonlinear, although gauge invariant; linear constitutive equations are indeed incompatible 
with Maxwell's equations and Galilean covariance. Further, we extend the Galilean sym- 
metry to the minimally-coupled Schrodinger-Maxwell theory, and to the coupled systems 
of Maxwell and nonlinear Schrodinger equations as proposed by Doebner and Goldin. 
Thus one has a consistent and fully Galilean covariant (but nonlinear) quantum electro- 
dynamics. 

There is a commonly held "folk belief" among physicists that Maxwell's equations, 
being Lorentz invariant, are not consistent with Galilean symmetry. However it is known 
(but not widely appreciated) that Galilean covariance, like Lorentz covariance, is a prop- 
erty of all four of Maxwell's equations for media in classical electrodynamics. More than 
25 years ago, Le Bellac and Levy-Leblond emphasized this in noting that the clash 
between Maxwell's equations and Galilean relativity occurs only in the constitutive equa- 
tions. Keeping the standard constitutive equations for the vacuum, and introducing some 
additional conditions for the electromagnetic fields, these authors derived two distinct 
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Galilean limits of Maxwell's equations: one in which Faraday's law is lost (called the 
electric limit); and another in which the displacement current is zero, that violates the 
continuity equation for charge and current densities (called the magnetic limit). Brown 
and Holland, in a recent discussion of these results, observed as follows: "It is noted that 
no fully Galilean-covariant theory of a coupled Schrodinger-Maxwell system (where the 
density and current of the Schrodinger field act as source of the nonrelativistic Maxwell 
field) is possible." § 

Dyson, in discussing an unpublished 1948 "proof" of Maxwell's equations by Feynman, 
remarks, "The proof begins with assumptions invariant under Galilean transformations 
and ends with equations invariant under Lorentz transformations. How could this have 
happened?" ||. Dyson's paper immediately provoked a heated discussion in the litera- 
ture; see, for example @, |5|, |6|. In particular, Vaidya and Farina ask in the title of 
their paper, "Can Galilean mechanics and full Maxwell equations coexist peacefully?" 
and answer forcefully, "No, they cannot." 

The fact that Maxwell's equations for media are both Lorentz and Galilei invariant 
sheds light on the mystery behind Feynman's insight. The conclusions of Brown and 
Holland, and Vaidya and Farina, depend on the implicit assumption of linear constitutive 
equations. In this letter we take a completely different point of view. We show how the 
choice between a Lorentz or Galilei invariant theory can be made by introducing a special 
class of nonlinearities in the constitutive equations, without modifying any of Maxwell's 
equations (or the continuity equation). Unlike Le Bellac and Levy-Leblond we are not 
forced by constitutive equations to choose between a privileged frame of reference (the 
ether) and Einsteinian relativity. 

Our interest in nonlinearity and Maxwell's equations is heightened by recent ex- 
perimental results in quantum optics, including optical squeezing and quantum non- 
demolition measurements, that suggest the fundamental importance of nonlinear con- 
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stitutive equations in media. In particular, we are motivated by recent experimental 
results demonstrating extremely slow light speed (as slow as 17 m/sec) in laser-dressed 
ultra-cold atomic media p] to look again at the question of Galilean-invariant electrody- 
namics in the context of nonlinear, coupled Schrodinger-Maxwell theories. The result is 
the framework for Galilean quantum electrodynamics described here. We would like to 
raise the question whether the constitutive equations in our Galilean class could be a good 
mathematical model for some real media, in some range of velocities and field strengths. 



2 Symmetry of Maxwell's equations for media 

A systematic investigation of symmetry of Maxwell's equations for media was reported 
in Refs. [f|, ||, where it was shown that these equations are covariant under the inho- 
mogeneous group of general linear transformations GL(A, R), that includes both Lorentz 
and Galilei transformations. General constitutive equations restricting this symmetry to 
Poincare symmetry were found. 

It is important for further analysis to use a system of units such as the SI system, 
that avoids incorporating c into the definition of the fundamental fields. Without such a 
choice it is not possible to see how nonrelativistic theory arises from a relativistic theory 
in the limit c — > oo. Maxwell's equations for media written in SI units have the form 
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VxE = -w 

V-B = 0; 

VxH= — +j, (2.1) 
V • D = p, 

where E is the electric field, D is the electric displacement, B is the magnetic induction, 
and H is the magnetic field; p and j are charge and current densities. The physically 



detectable fields are E and B, via the Lorentz force on a charged particle F = g(E+ v x B) . 
The fields H and D may be regarded as constructs used to describe (via the constitutive 
equations) how the directly observable fields are produced by charges and currents. 

As is well known, Maxwell's equations are Lorentz covariant; in particular, they are 
covariant under the space-time transformations |T^, [II], |T2| 

ar[l= 7 (x-v% x' ± = x ± , t' = 7 (t-^), (2.2) 

1 

y[\ -V 2 /C 2 ' 

with the corresponding field transformations 

B[=B h B' ± = 7 (B - iy x E) ± , 
E[ =E h E' ± = 7 (E + v x B) ± , 

H[ = H h H' ± = t(H- v x D)_l, (2.3) 
D\=D h D' ± = 7 (D + \v x H) ± ; 

II c z 

and the transformations for current and charge densities: 

v • j 

iy = 7(j - vp)||, j±=j±, P =l{p (2-4) 
It is appropriate to list here the Lorentz field invariants: 

h = B 2 -E 2 , J 2 = B ■ E; 

c 2 

h = B 2 - ^H 2 , h = H D; (2.5) 

h = B H - E D, J 6 = B D + H. 

c 2 

One can also verify that Maxwell's equations (2.1) are covariant with respect to the 
Galilean transformations 



E' = E + vxB, B' = B; 

H' = H-vxD, D' = D; (2.6) 
j' = j - P v > ft = P> 

which arise as the c — > oo limit of the Lorentz transformations (2.2) - (2.4). Both Galilean 
and Lorentz symmetries belong to the class of general linear transformations GL(4, R) 
admitted by Maxwell's equations (2.1). This, in turn, is a consequence of the incomplete- 
ness of the system (2.1): there are 8 equations for 12 unknown functions. The system 
must be completed by the constitutive equations, which are functional relations between 
vectors D, E, B, and H. 

A particular choice of the constitutive equations can reduce the symmetry of the 
system (2.1) to Lorentz or to Galilei. General constitutive equations that restrict the 
GL(4, R) symmetry of Maxwell's equations to Lorentz (Poincare) are reported in Refs. 
H In SI units, they have the form 

D = MB + ^NE, H = NB- ME, (2.7) 

where M and iV are arbitrary scalar functions of the Lorentz invariants Ji, I2 given in 
(2.5). Equivalently we may write 

B = RD + \QH, E = QD - RH, (2.8) 

where Q and R are arbitrary functions of the Lorentz invariants ^3, I4 from (2.5). 

Let us note that some well-known nonlinear theories like Born-Infeld electrodynamics, 
or Euler-Kockel electrodynamics that takes into account quantum-mechanical nonlinear 
effects (see, for example, |10| and references therein), correspond to particular choices of 



M and N in (2.7). 

The constitutive equations that reduce the GL(4,R) symmetry group of (2.1) to the 
Galilei group are 

D = MB, H = iVB — ME, (2.9) 



or, equivalently 

B = RD, E = QD — RH, (2.10) 

where M and N, Q and R are arbitrary functions of Galilean invariants. To demonstrate 
this, we first use (2.6) to see that the general constitutive equations D = f(B,E), H = 
g(B, E) must take the form D = MB, H = NB + N X E, where M, N, Aq are some 
scalar functions of Galilean invariants. Then substitution of the latter equation into (2.6) 
results in the condition Ni = —M. Finally, one shows that the Galilean field invariants 
are: 

h = B 2 , h = B - E; 

J 3 = D 2 , / 4 = H D; (2.11) 

h = B H - E D, J 6 = B D. 

One can, if one wishes, restrict oneself to constitutive equations in explicit (i.e., non- 
implicit) form, so that M and N depend on I±, I 2 , and Q and R depend on J 3 , J 4 . 

Note that the constitutive equations (2.9), (2.10) and the Galilean field invariants 
(2.11) are respectively the formal limits as c — > oo of the Lorentz invariant relations (2.7), 
(2.8), and (2.5). 

One obtains the standard Maxwell equations for the vacuum (where c~ 2 = eo/xo) by 
choosing N = l//i , M = in (2.7). But in the case of Galilean constitutive equations, 
letting M be a constant is not compatible with a nonzero charge density p. Therefore, we 
conclude that a consistent Galilean electrodynamics is essentially nonlinear. 

It is worth remarking for the case of relativistic constitutive equations (2.7), that 
letting M = A = constant, so that 

D = AB + 4^E, H = iVB — AE, 

c 2 

results in equations for the observable fields B and E that are independent of the magni- 
tude of A. In other words, constant terms A in the above constitutive equations have no 
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effect on observable fields, and therefore, without loss of generality one can set A = 0. In 
general, adding a constant to M in (2.7) has no effect on the observable fields E and B. 
This is also true for M in (2.9), which is why nontrivial Galilean constitutive equations 
are essentially nonlinear. 

The (linear) Galilean electrodynamics discussed in the literature is associated with 



the pre-Maxwell electrodynamics, with no displacement current [[[], |2|, [10], 112] , it is also 
called the "magnetic limit" of Maxwell's equations, and is the only one of the two limits 
discussed in [1] that is compatible with Schrodinger quantum mechanics. These equations 
in SI units have the form: 

<9B 

V-B = 0; 

VxB = /i j, (2.12) 
V • E = e ~V 

Let us note, that the approach used in to derive (2.12) from (2.1) actually embodies 
contradictory assumptions, moving back and forth between the two symmetries. First, the 
constitutive equations for the vacuum D = e E, B = /i H are imposed, immediately 
breaking the Galilean symmetry of Maxwell's equations (2.1). Then additional, non- 
invariant conditions about relative field strengths are assumed to justify the limit 1/c 2 — > 
in Maxwell's equations, and remove the displacement current. At the same time both eo 
and fio are kept non-zero, even though their product 1/c 2 has been taken to zero. Since 
there is no displacement current in (2.12), the continuity equation does not hold. The 
Galilean transformations for E and B in (2.12) are as in (2.6), while the corresponding 
transformations for p and j are 

j'=j, p = p - e /i v • j. (2.13) 
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In order to explore the consistency of the Galilean electrodynamics (2.1), (2.9), we 
have investigated two situations: first, the case M = 0, N = E • B (which requires p = 0); 
secondly, M — E ■ B, N — 1/pio- We were able to obtain some particular solutions, and 
this is a subject of our ongoing research. Let us note here that using the Galilei invariance 
of the system, one can construct (new) traveling wave solutions from (old) solutions by 
means of the following formulas: 

Enew(x) = Vold(x') - V X B oid (x'), B new (x) = B old (x')] 

H-new(x) = U old (x') + v x T) old (x'), T) new (x) = T> old (x'); (2.14) 

3new{x) = iold{x') + Vp old (x'), p n ew{x) = Pold{x'), 

where x' = (x', t') — (x — vt, t). 



3 A framework for a Galilei covariant quantum elec- 
trodynamics 

As is well known, any solution for vectors E and B of Maxwell's equations (2.1), can be 
represented in terms of potentials ($, A) 

B = V x A, 

E = -£-V». (3.1) 

for some $ and A. The choice of $ and A is not unique: new potentials A' = A + 
VA, $' = $ — dA/dt, containing an arbitrary function A, result in the same E and B 
(gauge invariance). The standard procedure of introducing electromagnetic interaction in 
quantum mechanics is minimal coupling, which is consistent with Galilean covariance of 
the free field equations. One obtains the coupled Schrodinger-Maxwell equation 

ik^- = — (-jW-eA)V + e^, (3.2) 
9 



where the electromagnetic fields E, B are obtained from $, A via (3.1), and are governed 
by (2.1). Gauge invariant current and charge densities entering (2.1) are given by 

p = W, J gi = ^Wip-NM] -— pA. (3.3) 
2im m 

Equation (3.2) is Galilean invariant, with corresponding transformations given in (2.6), 
and the transformations for vector-potential ($, A) given by 

A' = A, $' = $ - v • A. (3.4) 

The resulting fully Galilean covariant quantum electrodynamics is embodied in equations 
(2.1), (2.9) (with some concrete choices of M and N), (3.1), (3.2), (3.3). 

It is natural to ask about the coupling of nonlinear Maxwell theory with nonlinear 
Galilean invariant Schrodinger mechanics. A possible framework for such a quantum 
mechanics is given by the family of nonlinear equations proposed by Doebner and Goldin 
13| . Nonlinear terms of the form 



ihDAp V-j Ap j 2 j.Vp (Vpf 

—z ip + hD [ci hc 2 hc 3 — + c 4 — — + c 5 — — \ip (3.5) 

2 p p p p 2 p 2 p 2 

are added to the right side of Eq.(3.2), where j = (l/2i)['0V'0 — (Vfyip], and D,D' are 
diffusion coefficients. The Galilean-invariant subfamily of quantum theories is defined by 
C\ + c 4 = c 3 = 0. We are proposing consideration of the coupled system of Doebner-Goldin 
equations and nonlinear Galilean invariant Maxwell equations, as a consistent Galilean 
quantum electrodynamics. In this system, the quantities p, E, B, and 

jgi = jLupvip - - DVp- — pA. (3.6) 

2im m 

are gauge invariant. Eq.(3.6) gives the current to be entered in Maxwell's equations. Doeb- 
ner and Goldin further generalize the notion of gauge transformation to allow nonlinear 
gauge transformations, writing more general formulas for the gauge invariants p, J gl , E, 
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and B |L4]]. Recently Galilean invariant Doebner-Goldin equations were applied to de- 



scribe the dynamics of matrix D-branes [16 



In conclusion, let us remark that one can extend our approach further to include in 
the theory Galilei invariant spinor equations (see, for example, H). Another possible 
generalization of the quantum electrodynamics described here is to non-Abelian Yang- 
Mills fields. 
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